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SUMMARY 


The  turbulent  boundary-layer  drag  is  studied  for  production-type  or 
non-ideal  aircraft  sur faces.  Stemming  from  the  results  of  flight-test  data 
for  full-scale  late-model  aircraft;  it  has  previously  been  concluded  that 
there  are  no  apparent  Reynolds  number  effects  on  the  zero-lift  drag  at  flight 
speeds  below  the  drag-rise  Mach  number.  This  independence  of  boundary  layer 
drag  on  Reynolds  number  is  made  rational  by  reviewing  the  internal  mechanism 
of  the  incompressible  turbulent-boundary  layer  together  with  the  basic  em¬ 
pirical  results;  all  as  affected  by  surface  roughness. 

It  is  concluded  from  fundamental  concepts  that  the  effects  of  roughness 
elements  protruding  from  the  laminar  sub-layer  will  increase  the  vorticity 
and  thus  increase  the  momentum  losses  in  the  turbulent  boundary  layer.  This 
additional  momentum  loss,  or  drag,  due  to  roughness  has  been  shown  from  the 
empirical  evidence  to  offset  the  decreasing  skin  friction  drag  as  the  Reynolds 
number  increases. 

Correspondingly,  the  above  reasoning  is  extended  to  include  the  effects 
of  compressibility  in  the  turbulent  boundary  layer  up  to  where  the  roughness 
elements  themselves  become  locally  affected  and  no  longer  are  immersed  in  an 
incompressible  field.  As  a  consequence,  for  certain  conditions,  the  compres¬ 
sible  turbulent  boundary  layer  drag  for  rough  surfaces  may  also  be  considered 
independent  of  Reynolds  number. 

A  recipe  is  presented  for  estimating  the  boundary  layer  drag  contributions 
of  practical  aircraft  insulated  surfaces:  Given  a  surface  finish, an  equivalent 
sand-grain  roughness  is  determined  from  Table  I.  The  incompressible  boundary 
layer  drag  is  then  determined  from  Fig.  5  and  corrected  for  Mach  number  by 


HM-1129 


-2- 


Fig.  4.  The  resulting  drag  coefficients,  are  independent  of  increasing 

r 

Reynolds  numbers,  R  ,  for  R  greater  than  critical,  i.e.,  the  R  at  which 
equals  the  drag  due  to  pure  skin  friction,  C^.. 


0.472 


(13,14) 


At  R  less  than  critical  the  boundary  drag  is  that  of  pure  skin  friction  and 
e£ 

given  by  the  above  relation.  For  a  surface  of  given  length,  £,  measured 

streamwise  having  an  external  finish  equivalent  to  a  sand  grain  roughness, 

k  ,  such  that  k  /&  =  0.746x10”'’;  the  turbulent  boundary-layer  drag  estimate 
s  s 

shown  in  Fig.  3  may  be  prepared. 
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NOMENCLATURE 


Symbol  Definition 

x,  y  coordinates  parallel  and  normal  to  the 

surface  respectively 

u  local  velocity  in  the  x-direction 


M 


Mach  number 


tr 

p 

& 

K 

k 


local  boundary  layer  shear 
local  density 

turbulence  mixing  length 

mixing  length  constant 

mean  surface-roughness  height 


k 


s 


sand  grain  diameter 


Cf,  Cf 


local  and  average  skin  friction 
coefficient  respectively 


c 


f  5 


r 


local  and  average  rough-surface 
boundary  layer  drag  coefficient 
respectively 


A,B,  C,n 
T 


dT 


constants 

local  boundary  layer  temperature 

insulated  surface  or  boundary  layer 
condition 


Subscripts 


Reference 


w 


wall  or  surface 


00 


free  stream 


THE  TURBULENT  BOUNDARY-LAYER  VELOCITY  PROFILE-GENERAL  DISCUSSION 


The  purpose  of  this  section  is  to  extrapolate  the  known  effects 
of  surface  roughness  in  an  incompressible  turbulent  boundary-layer  to 
those  for  compressible  flow  regions.  The  mathematical  difficulties  en¬ 
countered  in  the  turbulent  boundary-layer  study  field  are  responsible  for 
a  variety  of  approximate  analyses  which  have  been  formed  around  two  or 
three  sets  of  original  data.  Although  semi-empirical  methods  abound 
for  describing  the  boundary-layer  phenomenon,  in  the  main  they  all  rest 
upon  a  set  of  fundamental  concepts  intended  to  ease  the  mathematical 
treatment  of  an  otherwise  unsolvable  problem.  It  is  important  to  realize 
what  the  limitations  are  of  these  fundamental  simplifications  in  order 
to  assess  the  results  predicted  from  their  use. 

Accordingly,  the  assumptions  and  accompanying  analysis  of  the 
internal  mechanism  of  the  turbulent  boundary-layer  are  reviewed.  One 
important  conclusion  made  is  that  due  to  the  emission  of  the  molecular 
viscous  effects  which  are  actually  present  at  or  near  the  wall  in 
boundary-layer  flow,  the  predicted  magnitude  of  local  velocities  occurring 
in  this  vicinity  are  expected  to  be  in  excess  of  those  physically  realized. 

Physically,  the  turbulent  boundary-layer  is  divided  into  three 
regions  of  differing  "viscous  action"  c.f.  Fig.  1. 

These  regions  are: 

a.  the  laminar  sub-layer  (molecular  diffusive  forces 
only) 

b.  the  buffer  layer  (molecular  and  eddy,  or  "fluid 
chunks",  diffusive  forces  of  comparable  influence) 

c.  the  turbulent  domain  (eddy  diffusive  forces 
predominate) 
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velocity  profile 


c.  eddy-visccus  forces 
predominate 

b. Jmolecular-viscPus  forces 
(plus  eddy- v±sc«ua  forces 

a.  molecular-visc pus  forces 


Fig.  1  THE  TURBULENT  BOUNDARY— LAYER  REGIONS  OF  VISCOUS  ACTION 


In  order  that  the  equations  of  motion  governing  the  momentum 
diffusion  or  energy  charge  through  the  incompressible  boundary  layer  be 
in  a  form  suitable  for  solution,  it  becomes  necessary  to  disregard  the 
effects  of  molecular  viscosity.  Therefore,  in  Fig.  1,  only  the  eddy 
forces  or  pure  turbulence  remain  to  affect  the  predicted  conditions  in 
the  boundary  layer. 

The  diffusion  of  the  eddies  create  energy  changes;  this  is 
due  to  the  momentum  changes  occurring  up  through  the  boundary  layer. 

The  result  is  a  local  resistance  force  or  apparent  shear.  The  individual 
eddies  involved  in  the  diffusion,  or  random  migrations,  are  considered  to 
be  solid  bodies;  therefore  at  each  point  in  the  boundary  layer,  the  resis¬ 
tance  force  is  then  proportional  to  the  square  of  the  local  velocity. 
Expressing  the  velocity  by  the  gradient  times  a  length,  the  local  apparent 
shear  is  described  by 


r 


dy 


du 

dy 


(l) 


where  Z  is  called  the  mixing  length  and  represents  the  turbulence  strength 
Prandtl  concludes  a  simple  form  for  the  mixing  length 
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as 


£  “  K-jjr 


(2) 


where  ■  constant.  Eq.  2  is  based  on  the  physical  observation  that  the 
vorticity  strength  decreases  down  through  the  boundary  layer  to  zero  at 
the  wall. 


Von  Harman 1 s  dimensional  analysis  of  the  mixing  length  concept 
results  in 


£  -  K„ 


du 

A 2_ 

,2 

d  u 

TT 


(3) 


where  “  constant.  Eq.  3  results  from  considering  the  eddy  migrations, 
at  each  point  in  the  boundary  layer,  normal  to  the  direction  of  the  mean 
flow  to  be  similar  when  the  sizes  of  the  fundamental  dimensions  of  length 
and  time  are  suitably  changed. 

It  is  repeated  that  both  Eqs.  2  and  3  result  from  disregarding 
the  energy  dissipated  by  molecular  viscosity. 

Prandtl  points  out  that  according  to  Von  Harman's  similarity 
theorem,  the  local  average  velocity  near  the  wall  must  be  of  the  form 


u  -  A(y+B)n  +  C  (4) 

with  A,  B,  C  and  n  are  constant.  Thus  from  Eq.  3 

£  =  K3  (y+B)  (5) 


where  y  is  the  center  of  similarity  at  each  point  in  the  boundary  layer, 
and  B  is  the  constant  which  shifts  the  reference  axis  on  which  £  m  0. 

When  Eq.  1  is  integrated  with  Prandtl' s  mixing  length  relation, 
Eq.  2,  considering  constant  shear  throughout  the  region,  there  results 
for  the  local  velocity 


HM-1129 


-7- 


u  *  J—  —  log  y  +  const.  (6) 

P 

Eq.  6  is  known  as  the  logarithmic  law  for  the  velocity  distribution 
throughout  the  turbulent  boundary-layer. 

In  Eq.  6,  when  y  =  0,  u  »  —  oo ;  this  is  contrary  to  the 
physical  condition  set  by  Eq.  2.  Therefore,  from  the  concept  of  similarity, 
B  >  0  in  Eq.  5j  and  the  reference  axis  of  similarity  lies  under  the  sur¬ 
face  of  the  wall.  At  the  wall  £ >  0  and  therefore  u  >  0.  Thus 
the  velocity  profile  represented  by  Eq.  6  is  seen  to  be  that  shown  in 
Fig.  2. 


Physical  case 
with  molecular 
viscosity 

Fig.  2.  TURBULENT  BOUNDARY  LAYER  VELOCITY  PROFILE 

Actually  it  is  known  that  the  molecular  viscous  forces  in  the  laminar 
sub-layer  will  halt  the  average  fluid  motion  at  the  wall;  therefore  it 
is  expected  that  all  turbulent- flow  velocity  profiles  at  or  near  the  wall 
will  have  actual  values  less  than  that  predicted  by  present  day  analyses. 
This  has  been  substantiated  in  all  measurements  that  have  been  compared  to 
analytical  methods  bassd  on  the  mixing  length  concept. 


THE  TURBULENT  BOUNDARY  OVER  A  ROUGH  SURFACE-INCOMPRESSIBLE  FLOW 


All  turbulent  boundary-layers  include  a  relatively  small 
laminar  sub-layer  adjacent  to  the  solid  wall.  This  region  behaves  as 
a  blanket  covering  such  surface  defects  as  crevices,  protruberances,  and 
roughness.  So  long  as  surface  irregularities  are  immersed  within  the 
low-velocity  flow  of  the  sub-layer,  the  over  all  effects  are  that  of  a 
smooth— wall  creating  pure  skin  friction  drag. 

When  surface  irregularities,  particularly  particles  of  rough¬ 
ness,  extend  outside  the  laminar  sub-layer  into  the  field  of  higher  mean 
velocities,  further  momentum  losses  occur  in  the  boundary  layer.  This 
additional  roughness  drag  together  with  the  pure  skin  friction  or  smooth- 
wall  effect  constitutes  the  total  boundary-layer  drag  experienced  by  most 
practically-constructed  and  finished  aircraft  components. 

To  date,  the  only  systematized  experimental  work  on  the  effect 
of  roughness  exists  for  incompressible  flow.  Nikuradse,  Ref.  1,  worked 
with  pipes  that  were  artificially  roughened  with  sand  grains.  His  experiment 
showed  that  for  a  given  roughness  there  exists  a  Reynolds  number  at  which 
the  measured  boundary  layer  drag  rises  in  value  above  a  trend  character¬ 
izing  a  pure  skin  friction  effect.  Furthermore,  with  increasing  Reynolds 
number,  the  deviation  from  pure  skin  friction  increases  until  a  constant 
total  drag  value  is  reached.  The  total  boundary-layer  drag  then  becomes 
independent  of  Reynolds  number  and  proportional  to  the  square  of  the 
velocity.  When  the  boundary-layer  drag  becomes  independent  of  Reynolds 
number,  fully  developed  roughness  flow  is  said  to  exist. 

The  increasing  roughness-drag  phenomenon  may  be  explained  in 
two  ways,  both  due  to  the  thinning  of  the  laminar  sub-layer  with  increasing 
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Reynolds  number.  The  first  is  that  local  flow  separations  are  induced 
as  portions  of  the  surface  emerge  into  the  region  of  turbulence;  this 
creates  a  pressure  drag  in  the  boundary  layer.  Apparently,  for  fully 
developed  roughness  flow,  the  additional  drag  increase  negates  the  de¬ 
creasing  effect  Reynolds  number  has  on  pure  skin  friction. 

Another  explanation  of  roughness  drag  is  to  consider  that  the 
particles  emerge  from  the  laminar  sub-layer  and  shed  additional  vortices. 
The  turbulence  intensity  is  increased,  but  the  mechanism  of  dispersion 
is  the  same  as  for  a  smooth e-wall.  Although  the  form  of  the  mixing 
length  equation  would  remain  unaltered,  the  reference  axis  of  similarity, 
B  in  Eq.  5>  would  change.  Correspondingly,  the  velocity  distribution 
for  roughness  flow  would  be  similar  in  shape  to  the  smooth  — wall  counter¬ 
part  but  of  smaller  local  values. 

Prandtl  and  Schlichting,  Ref.  2,  converted  Nikuradse's  results 
for  sand  roughened  pipes  to  flat  plate  flow.  The  boundary  layer  drag 
coefficients  for  fully  developed  roughness  flow  over  a  flat  plate  may 
be  represented  by  their  interpolation  formulas: 

,  -  2.5 

cf  -  (2.87  +  1.58  log10  £-  )  (7) 

r  s 

,  -  2.5 

Cf  -  (1.89  +  1.62  log10  f-  )  (8) 

r  s 

where  c^  and  are  the  local  and  average  coefficients  respectively, 
r  r 

The  velocity  distribution  in  the  turbulent  boundary  layer  with 
fully  developed  roughness  flow  was  derived  by  Prandtl,  Ref.  3,  from 
Nikuradse's  data,  and  is 
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-  8.48  +  5.75  log10  (9) 

s 

where  T'w  is  the  local  shear  at  the  wall  and  k  is  the  average  sand 

s 

grain  diameter.  Eq.  9  is  of  the  form  of  the  logarithmic  law,  Eq.  6; 
this  tends  to  substantiate  the  argument  for  the  increased  vorticity 
concept  in  explaining  the  mechanism  of  roughness  flow. 

THE  TURBULENT  BOUNDARY-LAYER  OVER  A  ROUGH  SURFACE-COMPRESSIBLE  FLOW 
Although  the  physical  mechanism  of  roughness  flow  can  be  ration¬ 
alized  from  two  independent  concepts,  viz.  induced  pressure  and  the  crea¬ 
tion  of  additional  eddies,  it  is  more  than  likely  that  both  properties 
contribute  to  the  roughness  drag.  By  virtue  of  the  fully  developed  rough¬ 
ness  drag  data  correlating  with  the  logarithmic  law  for  the  velocity  pro¬ 
file  it  is  concluded  that  the  increased  vorticity  effects  predominate. 
Furthermore,  the  scale  of  the  pressure  contribution  is  not  expected  to 
increase  unless  the  elements  of  roughness  are  subjected  to  critical  local 
flow  conditions  wherein  wave  drag  results.  Therefore,  in  extending  Eqs.  7 
and  8  to  boundary  layers  of  a  compressible  medium  and  supersonic  free 
stream  flow,  interest  is  first  centered  around  local  Mach  numbers  at  the 
roughness  elements. 

As  previously  discussed,  expressions  for  the  boundary  layer  velocity 
profiles  based  on  the  turbulence  mixing  length  concept,  to  the  exclusion  of 
molecular  viscosity  effects,  will  predict  larger  local  velocities  near  the 
wall  than  can  be  physically  realized.  Therefore,  any  of  the  existing  com- 
pressiblity  extensions  of  Eq.  1  can  be  employed  with  the  assurance  that 
the  predicted  effects  will  be  conservative  (larger  than  expected). 
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'  In  line  with  the  above  use  is  made  of  Van  Driest' s  compressibility 
extension  of  Prandtl's  apparent  shear  relation.  Ref,  7-  Van  Driest' s 
basic  assumption  is  that  Eq.  1  is  locally  valid  and  the  density  is  allowed 
to  vary  throughout  the  boundary  layer.  Proceeding  in  the  identical  manner 
it  can  be  shown  that  the  momentum  equation  in  the  region  of  roughness  is 


isin 


-1 


2A" 


es- 

oo 


B 


(B  + 


2372  *  Jsin 


-1 


B 


(B2  +  kk2)1/2 


where 


t-1 

2 

TW/T 
'  oo 


and  B 


i 

2  oo 


T¥/T 
'  oo 


1 


making  y  =  k  ,  u 
s 


(u  rough) 


max. 


and  Eq,  10  then  becomes 


jSxn  2k 


2  Urough 
u 


-  B 


(B  +  4A  ) 


C3  max, 

5375 


*  isin 


-1 


B 


(B^  +  4A  ) 


2072 


(11) 


The  energy  equation  is 


00 


jl«  _  1)  _H_  +tzl^  2  _u_  (1  _  _u_} 
T  VT  "7  u  2^<s>  u  u  ; 

CD  00  OO  00  00 


(12) 


For  present  airplane  considerations,  interest  in  the  boundary-layer 
compressibility  effects  does  not  extend  beyond  free  stream  Mach  numbers  of 
2.  Conservatively  assuming  that  local  Mach  numbers  below  0.4  at  the  tips  of 
the  roughness  elements  will  not  alter  the  incompressible  roughness  drag  effects. 
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it  is  possible  to  determine  the  first-order  maximum  allowable  degree  of 
roughness  before  compressibility  must  be  considered.  Inasmuch  as  the 
boundary-layer  temperature  influences  the  local  Mach  numbers,  two  extreme 
cases  of  cooling  and  heating  of  the  wall  are  considered;  viz.  temperature 
ratios 

^w/T  *  1.0  and  1.8  for  M  ■  2.0  . 

From  Eqs.  11  and  12,  the  local  boundary  layer  drag  coefficients  are 
tabulated. 


Mloc  " 

M  -  2.0 

00 

Tw/T 
'  00 

Uloc/u 

'  CO 

cf 

r 

1.0 

.213 

.0012 

1.8 

.260 

.0019 

It  can  be  shown  that  by  any  of  the  compressibility  factors  available  for 
ratioing  the  local  boundary  layer  drag  contribution  of  the  above  results 

to  the  incompressible  values,  and  by  substituting  the  equivalent  M  *  0 

c  1 

£  s  into  Eq.  7,  the  maximum  allowable  roughness  not  affected  by  com¬ 
pressibility  are  of  the  following  orders  of  magnitude: 


Mloc  '  -U 

M  -  2.0 

00 

Tw/T 
'  00 

^ s/£ 

1.0 

0(10“7) 

1.8 

0(10“5) 
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Therefore,  due  to  the  "smooth  —wall"  turbulence  diffusion  mechanism, 
the  same  compressibility  corrections  may  be  applied  to  rough-wall 
boundary-layers  for  free  stream  Mach  numbers  up  to  2  with  the  equivalent 
sand  roughness  not  exceeding  ratios  listed  in  the  above  table. 

Fig.  3. ,  shows  the  boundary  layer  drag  for  fully  developed 
roughness  flow  as  extended  by  the  additional  turbulence  concept.  All 
computations  were  carried  out  for  an  insulated  boundary-layer  and  equi¬ 
valent  sand  roughness  ratio 


ks/X  "  °‘76l^10~5 


The  dashed  curves  in  Fig.  3  represent  pure  or  smooth  —wall 
turbulent  skin  friction  according  to  Schlichtings 1  incompressible 
flow  relation. 


C 


f 


0.472 

(^lo  \>2'5 


(13) 


Q 

modified  for  compressibility  by  Fig.  4»  The  curve  of  vs  ^oo 

presented  in  Fig.  4  results  from  the  analysis  of  Rubesin  et  al.  Ref.  8. 
Included  with  the  compressibility  correction  curve  is  a  summary  of  some 
of  the  more  valid  experimental  results.  The  agreement  between  the  data 
and  semi-empirical  results  justifies  the  selection  of  this  particular 
curve.  According  to  the  results  of  the  analysis  in  Ref.  8  ,  there  is  no 
Reynolds  number  effect  for  the  skin  friction  parameters  combined  in  the 
above  manner. 

Analogous  to  Eq.  13,  the  average  drag  coefficient  for  fully 
developed  roughness  and  incompressible  boundary-layer  flow  described 
by  Eq.  8,  Fig.  5#  is  corrected  for  Mach  number  by  the  relation. 
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f 

M-0 


(1  +  M  2)0'467 
v  2  co  J 


,  M  <  4 

'  CD  ~ 


(14) 


of  Fig.  4. 


Therefore  within  the  bounds  of  the  present  discussion,  the 
boundary— layer  drag  coefficients  may  be  described  by  the  relations 


&  —2.5 

(2.87  +  1*58  log1Q  IT) 


(1.89  +  1=62  log1Q  ^-) 
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As  shown  in  Fig.  3,  it  is  assumed  that  the  turbulent  boundary- 
layer  flow  abruptly  changes  from  a  smooth — wall  condition  to  that  with 
fully  developed  roughness.  The  Reynolds  number  at  which  this  transition 
occurs  being  arbitrarily  defined  at  that  value  at  which  the  boundary-layer 
drag  coefficient,  ,  is  equal  to  the  smooth —wall  skin  friction  value,  . 


THE  DETERMINATION  OF  THE  EQUIVALENT  SAND  GRAIN  DIAMETER 


As  discussed  in  the  previous  section  the  application  of  the  method  for 
determining  the  boundary  layer  drag  coefficients  for  non-ideal  surfaces  is 
based  on  the  equivalent  surface  finish  measured  in  sand  grain  diameters. 

Table  I  provides  a  summary  of  the  types  of  surface  finish  correlated  with  the 
equivalent  sand  grain  sizes.  The  correlation  is  a  result  of  the  combined 
researches  of  Colebrook  and  White,  Hoerner,  and  Young;  Refs.  9,  10,  and  11, 
respectively. 


TABLE  I  EQUIVALENT  SURFACE  FINISH 


Type  of  Surface 

inches  x  leP 

kg 

inches  x  10^ 

Mirror  Finish 

0  -  .049 

0  -  .079 

Polished  Metal  or  Wood 

.05  -  .09 

.08  -  .159 

Natural  Sheet  Metal 

.1  -  .19 

.16  -  .319 

Optimum  Paint-Sprayed 

.2  -  .29 

.32  -  .479 

Average  Paint-Sprayed 

.30  -  .99 

.48  -  1.59 

Mass  Production  Paint-Sprayed 

1.0  -  10 

1.60  -  16. 

One  suggested  method  of  approach  is  to  determine  the  equivalent  sand-grain 

diameter  by  proceeding  backwards.  From  a  known  incompressible  of  a  given 

t. 

component  having  a  surface  finish  representative  of  future  designs,  the  kg  may 
be  determined.  With  but  a  few  samplings  of  each  component,  e.g.,  wing,  fuse¬ 
lage,  nacelles,  reasonable  estimates  may  then  be  made  for  projected  configurations. 
Such  contaminating  items  as  leaks,  isolated  clusters  of  rivet  heads,  and  the  like 
would  be  implicit  in  the  results. 


FIG.  4  TURBULENT  BOUNDS  LAVER  COMPRESSIBILITY  FACTOR 

INSULATED  CASE 


MACH  NUMBER- Mm 
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